function[ut] = get_ut(fday),

[iy, im, id, ih, min, sec] = datevec(fday+730486); % t is MJD2000

ut = ih+min/60+sec/3600;
